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Abstract: We explore the constraints imposed on higher curvature corrections of the 
Lovelock type due to causality restrictions in the boundary of asymptotically AdS space- 
time. In the framework of AdS/CFT, this is related to positivity of the energy constraints 
that arise in conformal collider physics. We present explicit analytic results that fully 
address these issues for cubic Lovelock gravity in arbitrary dimensions and give the for- 
mal analytic results that comprehend general Lovelock theory. The computations can be 
performed in two ways, both by considering a thermal setup in a black hole background 
and by studying the scattering of gravitons with a shock wave in AdS. We show that both 
computations coincide in Lovelock theory. The different helicities, as expected, provide the 
boundaries defining the region of allowed couplings. We generalize these results to arbitrary 
higher dimensions and discuss their consequences on the shear viscosity to energy density 
ratio of CFT plasmas, the possible existence of Boulware-Deser instabilities in Lovelock 
theory and the extent to which the AdS/CFT correspondence might be valid for arbitrary 
dimensions. 
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1. Introduction 

After several years of research, the AdS/CFT correspondence has become a standard tool 
to deal with non-perturvative phenomena in a plethora of systems that go far beyond the 
large N limit of J\f = 4 super Yang-Mills theory originally portrayed by Maldacena [1]. A 
recent application of this framework in the context of a gedanken conformal collider physics 
experiment was explored in [2] . A number of constraints for the central charges of 4d CFTs 
were shown to arise by demanding that the energy measured in calorimeters of a collider 
physics experiment be positive. For instance, every J\f = 1 SCFT must have central charges 
a and c lying within the window 1/2 < a/c < 3/2. It is well-known that gravity duals 
whose gravitational action is solely governed by the Einstein-Hilbert term lead to a/c = 1 
[3, 4]. In order to scrutinize CFTs with different values for their central charges, higher 
curvature corrections should be considered in the gravity side of the correspondence. 

The addition of a Gauss-Bonnet (GB) term for these purposes was explored recently. It 
was shown that the GB coupling should be correspondingly bounded. In the background 
of a black hole, the coefficient of this term, A, must lie within a given window in order 
to preserve causality at the boundary, i.e., if A lies outside the allowed region, boundary 
perturbations would propagate at superluminal velocities. This was originally shown in 5d 
where A < ^ [5, 6] and A > -7/36 [7, 8], the dual theories being 4d N = 1 SCFTs. It 
was later extended to 7d, —5/16 < A < 3/16 [9, 10, 11], the dual theories in this case being 
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6d (2,0) SCFTs. More recently, the result was extended to arbitrary dimensions with the 
result 1 [10, 11] 



(d-3)(3d- 1) 
4(d+l) 2 



< A < 



(d-4)(d-3) (d 2 -3d + 8) 
4 (d 2 - 5d + 10) 2 



where the actual existence of dual CFTs is less clear (see, however, [13]). Besides, quadratic 
curvature corrections were shown to arise in the string theory framework [14, 15, 16, 17]. 

Remarkably enough, both windows coming from such a different restrictions match 
exactly. This requires the use of the AdS/CFT dictionary to relate the 2- and 3-point 
functions of the stress-energy tensor (that account for the energy measured in the conformal 
collider setup) with the Gauss-Bonnet coupling. Such entry in the dictionary was originally 
obtained in the 5d case in [18], and recently extended to the 7d Gauss-Bonnet theory [9] (an 
alternative avenue was explored in [11]). These results seem to provide a deep connection 
between two central concepts such as causality and positivity of the energy in both sides 
of the AdS/CFT correspondence. Besides, these results provided an irrefutable evidence 
against the so-called KSS bound [19] for rj/s in quantum relativistic theories. This is due 
to the fact that the value for rj/s is corrected in presence of a Gauss-Bonnet correction to 
r)/s = j-(l — 4A). Since positive values of A are allowed, the shear viscosity to entropy 
density ratio, for such a SCFT, would be lower than the KSS value. 

As pointed out in [8], bounds resulting from causality constraints should not be a 
feature of thermal CFTs. The relation between causality and positivity must lie at a more 
fundamental level and, as such, should show up at zero temperature. Indeed, the upper 
bound on A in the 5d case was shown to come from causality requirements imposed on 
a scattering process involving a graviton and a shock wave [8]. It should be pointed out 
that the positive energy condition in CFTs used in the conformal collider setup is not self- 
evident at all. Hofman gave a field theoretic argument explaining why it holds in any UV 
complete quantum field theory [8]. He scrutinized deeper in the relation between causality 
and positive energy, and showed that there are indeed several bounds resulting from the 
different helicities both of the stress-energy tensor in the CFT side as well as of the metric 
perturbations in the gravity side. This was later extended to arbitrary higher dimensions 
in [10] (see also [11]). All these results provide a quite satisfactory understanding of the 
case of quadratic curvature corrections. 2 

A natural problem that immediately arises in this context is what happens in the 
case of higher curvature corrections? Whereas in the quadratic case any combination in 
the Lagrangian can be written as the Gauss-Bonnet term by means of field redefinitions, 
this is not the case for cubic or higher contributions. Indeed, for instance, there are 
two independent terms that can be written at third order [22]: the cubic Lovelock and a 
Riemman 3 combination. These two are very different in nature. While the cubic Lovelock 
term is the only third order curvature contribution leading to second order wave equations 
for the metric, in any space-time dimensions (d > 7), the Riemman 3 contribution leads 

The upper bound was found earlier in [12]. 
2 Related work has been pursued in [20, 21]. 
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to higher order equations of motion and, consequently, to ghosts. Furthermore, both 
terms lead to a different tensorial structure for the 3-graviton vertex, at least in flat space 
[22]. This means, a priori, that their inclusion should affect differently the corresponding 
tensorial contributions to the relevant 3-point functions of the stress-energy tensor. These 
tensorial structure is parameterized by two numbers, £2 and £4 [2] (that are related to 
the central charges of the CFTs). In the case of Gauss-Bonnet, all the results in the 
literature have £4 = 0, which is presumably related to the fact that the CFTs are SCFTs, 
in accordance with the intuition coming from the previous analysis in flat space. 

We will study in this paper the whole family of higher order Lovelock gravities. We will 
use the AdS/CFT framework to scrutinize these theories in arbitrary dimensions in regard 
of the possible occurence of causality violation. We shall address the more general case 
and explicitely work out the third order Lovelock gravity in arbitrary dimension. Results 
for higher order theories are contained in our formulas though some extra work is needed 
to extract them explicitely. The third order case has taught us that it is a subtle issue 
to determine which is the branch of solutions connected to Einstein-Hilbert and, thus, 
presumably stable (as far as we know, the discussion on Boulware-Deser instabilities [23] 
in Lovelock theory has not been fully established; we conjecture that there is at least one 
stable vacuum solution and, further, that it can be identified as the one connected to 
the Einstein-Hilbert theory). We also undertake the computation of gravitons colliding 
shock waves to seek for causality violation processes, and it seems clear to us that the 
identification of the relevant branch is of outmost importance. From the point of view of a 
black hole background, the relevant question is whether one can find an asymptotically AdS 
black hole with a well-defined horizon. This is a delicate problem that can be amusingly 
casted in terms of a purely algebraic setup. We show that the black hole computation and 
the shock wave one fully agree in Lovelock theory. 

We find the values of the quadratic and cubic Lovelock couplings that preserves causal- 
ity in the boundary. They define an unbounded (from below) region, whose intersection 
with the n = axis reproduces all the results previously obtained for Gauss-Bonnet. The 
region keeps its qualitative shape for arbitrary dimension but grows in width. As it happens 
in the quadratic case [10], one of the boundaries asymptotically approaches a curve that 
serves as the limit of the so-called excluded region where there are no well-defined black 
hole solutions. This is expected since that curve plays the role of an obstruction to the 
grow of the allowed region. The more striking behavior happens for the other boundary, 
which remains at finite distance when nothing seems to prevent its growing. Higher order 
Lovelock terms do not affect the value of r\js [24]. However, for positive cubic couplings, 
A can surpass the upper bound obtained in the pure Gauss-Bonnet case. This will push 
down the would be lower bound for r}/s. We will discuss these issues in the bulk of the 
paper. 

The structure of the paper is as follows. In Section 2 we present the main features of 
the Lovelock theory of gravity in arbitrary dimensions. We do this in the less familiar first 
order formalism in terms of tensorial forms, which is well-suited for most of the discussions. 
We analyze the singular locus of the theory and show that it divides the space of couplings 
in two well distinct regions. Section 3 is devoted to the discussion of black holes in Lovelock 
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theory. Even though there are some papers in the literature discussing this, we find that 
the analysis performed so far is uncomplete and pretty subtle. The theory has three 
branches of solutions, only one of which may be identified as the Einstein-Hilbert branch. 
The existence of a real solution in that branch displaying a black hole with a well-defined 
horizon is non-trivial. We formally solve the problem for arbitrary Lovelock theory. We 
focus in the case of third order in Section 4. There we use an amusing algebraic method to 
give a proper description of the different branches. Section 5 contains all the computations 
related to the possible violation of causality. We do this by two different ways. We actually 
prove that both coincide for any Lovelock theory in arbitrary higher dimensions. Section 
6 is devoted to the analysis of the results, to present our concluding remarks and to raise 
further questions. 

Note added: We were aware of an upcoming article on holographic causality in black hole 
backgrounds of Lovelock's theory by Jan de Boer, Manuela Kulaxizi and Andrei Parnachev 
[25]. Their article partially overlaps with some of our results. 

2. Lovelock theory of gravity 

Lovelock gravity is the most general second order gravity theory in higher dimensional 
spacetimes and it is free of ghosts when expanding about flat space [26, 27, 28]. This means 
that these theories possess the same degrees of freedom as the Einstein-Hilbert Lagrangian 
in each dimension. The action for this theory can be written in terms of differential forms 3 
as 

[— ] 

I= Y^kf [k) ( } 

^ d-2k J y 1 

where is the exterior product of k curvature 2-forms with the required number of 
vierbeins to construct a (i-form, 

n {k) = e fl ... fd R hh A • • • A Rhk-if2k A e hk+i-f d ? ( 2 .2) 

where eh—h is a short notation for &h A . . . A e^ k . For convenience, we will also use from 
now on the expression Rh-hk = Rhh a • • • A Rhk-ihk_ The ze roth and first term in (2.1) 
correspond, respectively, to the cosmological term and Einstein-Hilbert action. These are 
particular cases of the Lovelock theory. It is fairly easy to see that cq = -j? an d c\ = 1 
correspond to the usual normalization of these terms. 4 

3 We make use of the language of tensorial forms: out of the vierbein, e a , and spin connection, uj a b , 
1-forms, one can construct the Riemann curvature, R a b , and torsion, T a , 2-forms: 

R\ = du,\ + a," A U \ = ijiV dx» A dx v , T a = de a + c\ A e" . 

4 We set, without any lose, 167r(d — 3)! Gjv = 1 to simplify our expressions. In tensorial notation, thus, 
the cosmological constant has the customary negative value 2 A = — (d — l)(d — 2)/L 2 in d space-time 
dimensions. 
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If we consider this action in first order formalism, we have two equations of motion, 
for the connection 1-form and for the vierbein. If we vary the action with respect to the 
connection the resulting equation is proportional to the torsion 



k Cfc 



£ab = e abfs ...f d JZ^k ( Rh "' hk A e /2 *+ 1 -^" 1 ) a T^ d = . (2.3) 



k=l 



Thus, we can safely impose T a = as in the standard Einstein gravity. The second equation 
of motion is obtained by varying the action with respect to the vierbein, 



£a = eo/i-/rf_i Yl Ck ( Rh '" hk A ' hk + 1 ~ fd ) = • (2.4) 
fc=l 

If we further assume that ck is the highest non- vanishing coefficient, i.e., Ck>K = 0, this 
can be nicely rewritten as 

^•••/d-i Hi)" A • • • A ?m~ lhK A ^ hK+1 "- fd - 1 = , (2.5) 

where J-^ = R ab — Aj e ab , which makes manifest that, in principle, this theory admits K 
constant curvature ('vacuum 1 ) solutions, 

-pah = R ab _ A . e ab = Q (2 _ g) 

Plugging this expression into (2.4), we found that the K different cosmological constants 
are the solutions of the K-th order polynomial 

K K 

T[A]^^c fe A fc = c^[](A-A i ) = 0, (2.7) 

k=0 i=l 

each one corresponding to a ' vacuum 1 you can consider perturbations about. The theory 
will have degenerate behavior and possibly symmetry enhancement whenever two or more 
effective cosmological constants coincide. This is accounted for by means of the discriminant 

A = J](Ai - A,) 2 . (2.8) 

i<j 

Since we will later focus in the third order Lovelock case, let us discuss it in further detail. 
The lowest dimensionality where this term arises is 7d. Consider the following action, 



16irG N J V y 



(2.9) 



where the quadratic and cubic Lagrangians are 

£ GB = R 2 - iR^R^ + R^paR^ , (2.10) 

l2RR ttv R (lv + 24Ri* ua PR ail Rp u + WR^R^R^ 
+ 2AR^R aPvp R l f + 8R^ ap R a i a R pa ^ + 2R a p pa R^R^ v . (2.11) 



- 5 - 



It can be casted in the language of our previous discussion as (02 = XL 2 and C3 = ) 
/uL 4 f /„„w^ d-6 3A 



3(d-6) 

d-6 3 

+ 



d-4 n L 2 



R abcd A e e/ 



i? ab A e cde/ + 



d-6 3 



d — 2 [iL A d fiL e 

whose equations of motion can be written as (we set, again, T a = 0): 



(2.12) 




Figure 1: Singular locus for cubic Lovelock gravity. The shaded region is the domain where A > 0. 
At the (dotted) vertex, A = 1/3 and ii + (l/3) = /i~(l/3) = /x = 1/9, the theory becoming maximally 
degenerated (the symmetry enhances at that point to the full AdS group). The (dotted yellow) 
point, belonging to the locus, is nothing but the singular locus of Gauss-Bonnet gravity, A = 1/4. 
The point A = /i = naturally belongs to the curve: it corresponds to Einstein's gravity. 

( „... /(/: ...,,„ „ (i?° 6 - Ai e ab ) A [R cd - A 2 e cd ) A (i? e/ - A 3 e e/ ) A e 9 ^ 9 ^ = . (2.13) 

The values of the effective cosmological constants are complicated functions of the couplings 
that are not really important. They must satisfy, 



5> 



fiL 2 ' 



i<j 



i<j<k 



or, equivalently, 



^ (L 2 A) 3 + A (L 2 A) 2 + (L 2 A) + 1 = 
3 



(2.14) 



(2.15) 



There is always (i.e. for any A and //) at least one real cosmological constant. 5 The theory 
will have degenerate behavior whenever the discriminant (2.8), which can be nicely written 
in terms of the couplings of the action, 

1 



A = ~ (12 A 3 - 3 A 2 - 18 X/j. + n(? n + 4)) 



(2.16) 



J This is to be contrasted with the case of Gauss-Bonnet where 1 — 4 A should be positive for the theory 
to have an AdS vacuum. 
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vanishes. The singular locus is T : //(A) = /U + (A) U /i (A), with 

M ± (A) ^ A _2 ± 2 (1 _ 3A)3 / 2 (217) 

It is convenient, for later use, to write the discriminant as 

A = -3 ( M - /i+(A)) (/x - /i-(A)) , (2.18) 

since it makes clear that A < if fi > /i + (A) or fj, < /z _ (A), A = if \i = /i + (A) or 
[i = /i _ (A), and A > 0, elsewhere. We will see in what follows that the singular locus T 
plays an important role. Notice that it does not depend in the space-time dimensionality 
(for d>7). A similar analysis can be readily performed for higher-order Lovelock theories. 

3. Black holes in Lovelock theory 

Black hole solutions can be obtained by means of an ansatz of the form [29, 30, 31] 

ds 2 = -N% f(r) dt 2 + — + T1 dY? d _ %c , (3.1) 

where ^ 

dSd-V = 1 _ d f p 2 /L 2 + P 2dn d-z > ( 3 - 2 ) 

is the metric of a (d — 2)-dimensional manifold of constant curvature equal to (d — 2){d— 3)a 
(a = 0, ±1 parameterizing the different horizon topologies), and d^l 2 i _ 3 is the metric of the 
unit (d — 3)-sphere. For a = 0, which is the case we are going to be dealing with in this 
paper, using the natural frame, 

e° = N #v ^)dt , e i = __ * dr, e = jdx a , (3.3) 

v/( r ) L 

where o = 2, . . . , d — 1, the only non- vanishing components of the spin connection read 

lo° - f ' {r) e° oj 1 - e a (3 A) 

-i- 2 ^ W) e , « B - r e , (3.4) 

for torsionless space-time. The Riemann 2-form, R ah = dco ab + uj a c A co cb , 

BOX = * f » {r) e A e l R °a = _fp_ e A( ,a 

2 2r 



R la = _f(r) e l Ae a j ^afe = g a Ae fe _ ,35) 

2r 

(We omit in what follows, for the sake of clarity, the argument of the function.) If we insert 
these expressions into the equations of motion, we get 6 

£ = => Y, c k { -^(r(f k )' + (d-2k-l)f k )=0, (3.6) 

fc=0 



6 Recall that, for later use, we introduced K as the natural number, K < [^"H, that labels the highest 



non-vanishing coefficient, i.e., c^^k = 0. 
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that can be nicely rewritten as 

d 



dlogr 



+ (d-l) 








(3.7) 



where g = —f/(Ar 2 ) is basically the quotient between the black hole function and the 
vacuum solution (we are naturally assuming that the vacuum is given by one possible 
cosmological constant emerging from the analysis of the previous section). This can be 
easily solved as 

K 

K 



(3.8) 



fc=0 



where k is an integration constant somehow related to the mass of the spacetime 7 . It must 
be positive in order for the space-time to have a well defined horizon. To see this we can 
rewrite the precedent equation as 

X>AV = -£i-^, (3-9) 



k=i 



and realize that the equation admits a vanishing g only when r = r + = (kL 2 )^^ . We will 
thus reshuffle this relation and write everything in terms of r + . Notice that (3.8) leads 
to K different branches, each one associated to each of the cosmological constants (2.7). 
Just one branch has a horizon there. This can be seen since the polynomial root g = has 
multiplicity one at r = r+. In order to have higher multiplicity it would be necessary that 
also ci vanishes, but c± = 1 in our treatment (it corresponds to the Einstein-Hilbert term). 
For instance, in Gauss-Bonnet gravity there are two branches that read 



r 

2L 2 X 



1 ± 



N 



1 - 4A 1 



.d-l 



„d-l 



(3.10) 



Only one of the solutions, /_ , has a well defined horizon, the other one displaying a naked 
singularity at the origin. Remarkably, the solution with the well defined horizon is the one 
connected to the standard Einstein-Hilbert gravity, in the sense that it reduces to it when 
A -»• 0, 



f- 



2L 2 A 



1 



1 - 2A 1 



.d-l 



-d-l 




.d-l 



r.d-1 



This black hole has a well defined temperature 

T=N />+) _ V-1)N# r+ 
* 4vr 4vr L 2 

its entropy still satisfying the usual area law [32, 33, 34, 35] 

r d~2 



(3.11) 



(3.12) 



(3.13) 



7 For a 7^ equation (3.6) remains formally the same except for the shift / — > / — a; thus, an analogous 
redefinition of g, f = a — Ar 2 g, leads again to the relevant solution of (3.8). 
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where we have divided by the volume of the flat directions. The mass can be obtained as 



m = 



I 



Tds 



f 

Jo 



■r+ 



dr+ 



(d-2)N# f 



4-1 
+ 



(3.14) 



16vr 



L 2 



This is general for any Lovelock black hole of the type considered here, since the only 
term that contributes to the temperature is the Einstein-Hilbert one and the area law for 
the entropy has also been found to be valid in general for Lovelock black holes with flat 
horizons [36]. 

In general (3.8) has K solutions, one associated with each possible value of the effec- 
tive cosmological constant, and these are continuous functions since the roots of a poly- 
nomial equation depend continuously on its coefficients [37]. Rescaling the coefficients 
Cfc = afc L 2k ~ 2 and calling x = L 2 A g, 



where, of course, x = x{r). When r — > oo this is nothing but the equation for the (di- 
mensionless) cosmological constant (2.15). The different couplings afc>i fix the shape of 



Figure 2: The polynomial p[x; rj] for the case K = 2, i.e., cubic Lovelock theory (with A = 0.2 
and fx = 0.01), in arbitrary space-time dimension, for different values of the radius ranging from r+ 
to oo. The depicted points give x(rt) for the relevant branch. 

the polynomial p[x;r]. While varying r from r + to oo (see Figure 2), the curve translates 
rigidly along the vertical axis since the inhomogeneous coefficient ao = oo( r ) gives the point 
where these curves intersect it. It is a monotonically increasing function of r. Since p'[x; r] 
is positive at x = 0, the function x(r) is decreasing close to r + . But it should be monotonic 
since ao(r) is so and the remaining coefficients are frozen. Given that the cosmological 
constant is negative, g must be a monotonically increasing function of the radius. Hence, 




(3.15) 



PM 




X 
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from (3. 



r+9 



(d-1) r« 
L 2 A 7 C 



K 



-1 k-1 



k=l 



> , Vr G [r+,oo) 



(3.16) 



Then, the metric is regular everywhere except at r = and whenever the relevant solution 
coincides with any of the spurious ones since, in such case, 



K 

^2kc k A^g 

k=l 



k-1 



. 



(3.17) 



The value of r at which this happens exhibits a curvature singularity that prevents entering 
into a region where the metric (the function g) becomes complex. Also, when taking the 
limit r — > oo, (3.16) implies 



K 

^fccfcA^ 1 > . 

k=l 



(3.18) 



This vanishes when A = 0, since these are simultaneously zeroes of a polynomial, T[A], 
and its derivative, T'[A], and so higher multiplicity (at least second order) zeros of the 
original polynomial. However, the relevant branch can either be one of the degenerate ones 
or not and so not all the A = points will be relevant for our discussion. We shall refer to 
this later. Eq.(3.18) gives the derivative of the original polynomial at infinity. The slope of 
the polynomial must be non-negative for g to be monotonically increasing. This condition 
protects the relevant branch from Boulware-Deser-like instabilities [23]. 

We can further determine the asymptotic behaviour of g at infinity, g ~ \ — Rr d ^ l /r d ~ 1 . 
If we plug this expression in (3.16) we get 

-i 



1 

Ua 



■ K 

£ 

k=l 



kc k A K 



k-1 



(3.19) 



and, as long as the polynomial factor is positive and the cosmological constant negative, k is 
positive. The well-behaved black hole behaves asymptotically as the usual (positive mass) 
AdS black hole from Einstein-Hilbert gravity. The addition of higher curvature corrections 
does not change the qualitative behavior of the solution in any meaningful way. 



Im[x] 

0.2 r 



Re[x] 



Figure 3: The AdS black hole solution of the Einstein-Hilbert theory. 

This algebraic reasoning leads us to an alternative formulation: we can visualize each 
branch as a segment corresponding to the continuous running of the roots (in the complex 
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plane) while r goes from the boundary to the horizon. For instance, the solution of the 
standard Einstein-Hilbert gravity with negative cosmological constant is described by the 
root of the linear polynomial (see Figure 3), 



PEu[x;r] = 




(3.20) 



The Gauss-Bonnet case, instead, presents a richer structure. It is given by the roots of the 
quadratic polynomial 

PGB[x;r] = U-^ztJ +x + Xx 2 = . (3.21) 

The discriminant is A(r) = 1 — 4ao(r) A, which immediately suggests that there are two 
different regions. When A < 1/4, A(r) > Vr 6 [r+,oo), while for A > 1/4, there is a 

Im[x] 



oo r t 



Re[x] 



Figure 4: The AdS black hole solution of the Gauss-Bonnet theory for A = —0.4. There are two 
real branches. The green one is the branch connected to the Einstein-Hilbert AdS black hole, the 
one we named /_ in (3.10). 

value r+ £ [r + ,oo) where A(r*) = (see Figures 4 and 5). For r > r*, A(r) < and 
the solutions become complex (see Figure 6) behind the curvature singularity. The critical 
value A = 1/4, which is the singular locus of the Gauss-Bonnet theory, leads to r* = oo, 
the curvature singularity disappears but the theory has no linearized degrees of freedom 
around the AdS background. This quite simple behavior for different values of the GB 
coupling serves as a preliminary exercise to clarify the algebraic approach in higher order 
theories. The situation, indeed, gets more involved in the third order Lovelock theory. The 
relevant polynomial is 

PLU [x;r] = (^l- r -±-^ +x + Ax 2 + |:e 3 = , (3.22) 

whose discriminant is A(r) = ao(r) A (6/i — 4A 2 ) — 3a 2 , (r) /i 2 + A 2 — 4/i/3. As a function of 
the radius, A(r) = spans a family of curves in the (A, /x) plane that go from fi = | A 2 
(namely, A(r+) = 0) to the singular locus T depicted in Figure 1. It is not difficult to see 
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Im[x] 



Re[x] 



Figure 5: The AdS black hole solution of the Gauss-Bonnet theory for A = 0.2. For A > 0, the 
(still real) red branch jumps to the left of the green one, without flipping its orientation. For bigger 
values of A, it slides to the right, approaching the green branch. Needless to say, their infinities 
collide at A = 1/4. 



Imlxl 



Re[x] 



Figure 6: The AdS black hole solution of the Gauss-Bonnet theory for A = 0.26. For A > 1/4 the 
branches become complex conjugates to each other behind the curvature singularity (r = r*). The 
critical value, indeed, is the singular locus of the theory (see Figure 1). 



that for intermediate values of r, the curves look roughly like T, the singular vertex sliding 
up the parabola \i = A 2 (see Figure 7). Thus, the (A, y) plane has two regions, M.^\ where 
there is a value r* £ [r+, oo) such that A(r*) = 0, 



M {+) = i A < , n + (\) < fj, < ^A 2 j U jo < A < 



27 ' 4 



A 2 < n < y + (X) 



<27<A<-, At" (A) <y<y + (X) 



c - , //<//-(A)j> U <|A> — , y< -A 2 



(3.23) 
(3.24) 



and one region, M^ 2 \ where there are two values G [r + , oo) where A(r^) = 0, 



l ^ <X -l' I A2 <^<^(A)| U [a>±, ^A 2 <^<A 2 



(3.25) 
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Figure 7: Vanishing locus of A(r) for different values of r, spanning from oo (thickest red and 
blue curve in the left) to r = (20/17) 1 ^ d ~ 1 ' r + (thinest red and blue curve in the right). The vertex 
of the curve slides up along the parabola /i = A 2 (dashed curve). Points in the brown region (as 
the one depicted in yellow) belong to two curves. The black curve is fj, = |A 2 (that corresponds to 
A(r+) = 0), which bounds the regions colored in light blue, where A(r) only vanishes for one value 
of the radius between infinity and the horizon. The zoom in the right allows to better understand 
the structure where the different regions merge. 

Everywhere else, A(r) does not vanish for real values of r (unless they are hidden by the 
horizon). The previous analysis suggests that the space of parameters is divided into dif- 
ferent regions that should be treated separately. Depending on the sign of the discriminant 
we will have one or three real solutions (cosmological constants). As we shall see, for A > 
we have three real solutions while for A < we will have just one. 

The algebraic approach portrayed in this section can be extended to arbitrary higher 
order Lovelock theory in higher dimensions. It can be surely dealt with analytically up to 
the fourth order gravity while higher order cases may require numerical analysis since the 
quartic is the highest order polynomial equation that can be generically solved by radicals 
( Abel-Ruffini theorem) . We will focus from here on in the case of cubic Lovelock theory. 

4. Black holes in third order Lovelock theory 

As discussed in the preceding paragraphs, most of the information needed to clarify the 
existence of black hole solutions for different values of the Lovelock couplings relies in the 
behavior of a cubic polynomial. There is a very convenient way of parametrizing the cubic, 
suitable for analyzing the relevant branch of the black solution in terms of the geometry 
of the polynomial. We can characterize a cubic in terms of four fundamental parameters 
[38], 5, h, xn and (see Figure 8). Let us show how this works in our case. 

If we start from our cubic polynomial, Pll 3 [x] r] = ^x 3 + Ax 2 + x + ao( r )> a shift on 
the variable x = z + xn, where xn = — jj, leads to a polynomial, j>ll 3 [z; ?"], that is known 
as the reduced cubic. N is the point of inflexion, and it can be shown that h is a simple 
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Figure 8: Geometrical meaning of the parameters of the cubic. Recall also that the cubic is 
symmetric with respect to the point of inflexion. 



function of 5, namely 

h = --fi5 3 , where 5 2 = . (4.1) 

3 \i l 

Thus, the shape of the cubic is completely characterized by the parameter 5. Either the 
maximum and minimum are different (5 2 > 0), or they coincide at N (5 2 = 0), or there 
are no turning points (5 2 < 0). We choose the sign of 5 so that 8 > corresponds to the 
situation depicted in the figure (fj, < 0) and h (when S is real) is a positive quantity, 

5 = " * \/A^ => h = ^ (A 2 - M) 3/2 • (4.2) 

Consider the usual form of the cubic equation j?ll 3 [2;;^] = 0, with roots a, j3 and 7, and 
obtain the reduced form by the substitution x = + z. The equation has now the form, 

f *z 3 -»S 2 z + y N (r) = 0, (4.3) 

with roots a — xjy, (3 — xn and 7 — xn- The parameter yAr(r) is obviously the only one 
that depends on the radial coordinate through ao(r), 

2/7v(r) = a (r) + ^ (2A 2 - 3/x) . (4.4) 

This form allows us to use the identity 

(p + vf - 3p q(p + g) - (p 3 + q 3 ) = o . (4.5) 
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Thus z = p + q is a solution where 

pq = 5 2 and p 3 + g3 = _3yjV (46) 

Solving these equations by cubing the first and substituting into the second, and solving 
the resulting quadratic in p 3 gives 




The discriminant of the polynomial reads 

A(r) = -3/i 2 {y% - h 2 ) = -3fi 2 (y+ yJ) , (4.8) 

where y± = y^ ± h are the y-coordinates of the maximum and the minimum respectively. 
A(r) acquires a neat geometrical meaning in the light of this expression. We can see 
that the sign of the discriminant is determined by the position of the maximum and the 
minimum with respect to the x-axis. When h G C, y+, y_ £ C, this corresponding to the 
case /i > A 2 . The singular locus A = corresponds to y± = at r = oo (see Figure 1). 

We can treat separately the different regions. For y 2 N > h 2 (y+y- > 0), there is just 
one real root that can be easily obtained as 8 




For yjj = h 2 (y+y- = 0), there are three real roots two of them equal. The roots are 
a = xn + 25 and /3 = 7 = xat — <5, where the sign of 5 depends on the sign of y^ and has 
to be determined from 




(4.10) 



If UN = h = 0, then 5 = 0, in which case there are three equal roots at x = xj\[. 

For yjj < h 2 (y+y~ < 0), all three roots are real and distinct. The easiest way to 
proceed, without having to find the cubic root of a complex number, is to use trigonometry 
to solve the reduced form with the substitution z = 25 cos 9, that gives 

cos3fl=^. (4.11) 

8 The complex branches being, 
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The three roots are therefore given by 



a = xn + 25 cos 9 , 

P = x N + 25cos(6 + 2tt/3) , (4.12) 
ry = x N + 2(5cos(6» + 47r/3) . 

The important point to notice here is that, as y/y varies (or equivalently r since yAr(r) 
is monotonic), the angle can run from zero (where /3 = 7 and a is the "real" branch 
parameterized by (4.9)) to tt/3 (where = 7 and /3 is the "real" branch parameterized by 
(4.9)). We can follow the real root(s) as yjy changes. For y^ > h we have just (4.9) as a 
real root. For y/y = h we have a = xat + 25 and /3 = 7 = xtv — 5. For — /i < y/y < h, 
the angle is monotonic with yjy going from 6 = to tt/3. Beyond that point, two roots 
become complex again but not the same two that were imaginary for y/v > h. Thus, we 
have f3 = xn — 25 and = 7 = x^ + 5 for y^v = —h, the latter two roots becoming complex 
for y^v < —h where j3 corresponds to (4.9). Each of these solutions is continuous with y^ 
and the same can be checked with the other parameters of interest, [i and A. The only 
presumably singular point is /1 = where the degree of the equation changes, but we can 
take a well defined limit where one of the solutions diverges and the other two coincide 
with those of the quadratic polynomial. These solutions are continuous, but the 'always 
real branch' as usually parameterized in (4.9) is discontinuous since there is an interval of 
values of x that cannot be taken by the a nor the /3 branches (see Figure 14). From now 
on we will refer to as a (/?) the branch real for yjv — > 00 (—00). 

Just by analyzing the shape of the cubic we can understand which one is the branch 
connected to the horizon (x = for r = r + (ao = 0)). We have given values for 00 and oq, 
that is, the polynomial at the origin and its first derivative there. The first derivative of 
the polynomial at x = is always 1, and so the polynomial is growing at that point. Also, 
we know the sign of x^r when A and [i are given. We can distinguish several different cases. 
We know from our previous analysis that we can take a representative point in each of the 
relevant regions and analyze the behavior of the solutions. 

For [i > A 2 (5 2 < 0) there are no turning points (or they coincide at the inflexion point) 
and the polynomial is monotonically growing. There is just one real branch of solutions for 
all values of r and it has then no discontinuity. The same happens in the region contained 
in between the curve n = 3/4A 2 and fi = A 2 for A < and in between \x = /^ + (A) and 
fi = A 2 for A > 0. In this whole region we have A(r) < Vr £ [r + , 00). This happens for 
any point (A,/i) lying in the upper white part of Figure 7. Let us call this region 

M {0) = I A < , /i > ^A 2 } U |o<A<^, / u> M +(A)| U I A > ^ , /x > A 2 j . (4.13) 

In this case, the roots of (3.22) behave the same all along the radial flow; there is a single 
real branch which is the relevant one for asymptotically AdS black hole solutions with a 
well-defined horizon, and two complex conjugate unphysical branches (see Figure 9). The 
/x = A 2 line for A > 1/3 is excluded from M.^ since in this case the inflexion point is 
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Figure 9: A(r) < Vr 6 [r + ,oo) corresponding to Ad^K There is a single branch which is real. 
As we will see later, it is the relevant branch for the black hole solution. 



situated above the x-axis, 



Vn(oo) = -f\-- \ > I) 



(4.14) 



with xn < 0, and then, when the inflexion point crosses the x-axis (r = r*, where A(r*) = 
0), a (naked) curvature singularity shows up even if the spacetime is perfectly regular for 
r / r*. 

The lower white part of Figure 7, instead, has A(r) > Vr € [r+, oo). Let us call this 
region M. 



(0) 



(o) 



A < , ix-{\) <f x< } U <J < A < — , M -(A) < /x < -X 2 } . (4.15) 



In this case, again, the roots of (3.22) behave the same all along the radial flow (see Figure 
10); there are three real branches but only one is relevant for asymptotically AdS black 
hole solutions with a well-defined horizon (see Figure 10). This is reminiscent of the Gauss- 
Bonnet case where there are two real solutions for A < 1/4. One of the branches diverge 
as [i — > 0. 



• — — it 



-2 



/"+ OO 



Figure 10: A(r) > Vr € [r + ,oo) corresponding to . This situation is very similar to 
Gauss-Bonnet gravity. One of the solutions diverges as n — > 0. 
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For fx < A 2 we have in general two values of r (or y^) for which A(r) = 0, but they 
can be included or not in the interval r £ [r+,00) (in the previously analyzed regions 
these values lie outside this interval). For fi > (6 < 0) and A < (xn > 0), both the 
maximum and the minimum are located at positive values of x and so the branch connected 
to the horizon has no problems of reality or continuity even if eventually (in the subregion 
contained in Ai^) A(r*) = for some value r* G [r+,00) (see Figure 11). 



Figure 11: < and A(r + ) > corresponding to 

For n < (5 > 0) the only growing part of the polynomial is the one in between the 
minima and the maxima, corresponding to the branch 7. The maximum is located to the 
right of the origin x = and does not pose any problem regarding the relevant branch, 
but the minimum is located in negative values of x and then we must have ?/_(oo) < 
(un (co) < h) in order to have a real cosmological constant. The singular curve y_ = 
corresponds, in this region, to fi = /i~(A). The region with real cosmological constant 
corresponds to while that with complex cosmological constant corresponds to the 

region (see Figure 12). 



(1.6 



0.4 



0.2 



• < • 



Figure 12: Aoo < and A(r + ) > corresponding to A4^~K Regarding the behaviour of A(r) 
this region seems similar to the previous one but in this case one of the degenerate branches when 
A(r) = is the one connected to the horizon. Thus, there is a real cosmological constant but it 
does not correspond to the branch with horizon. 

The 7 branch for fj, < is continuously deformed into a or (3 when crossing the /i = 
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line. For A < the /3 branch diverges and 7 and a interchange their roles. For A > is a 
the diverging branch and 7 is deformed into f3. For fi = we can identify the remaining 
branches as the solutions for the Gauss-Bonnet case. 

The remaining case to be discussed is fi > (6 < 0) and A > (xn < 0). In this 
region (except as already discussed for fi > A 2 ) the two critical points are located at 
negative values of x, and they can have positive or negative y-values. For the subregion 

3 r 




-2 - 



-3 L 

Figure 13: > and A(r+) < corresponding to M^ + \ 



contained in see Figure 13. As in the previous case, in order to have a real value 

for the relevant cosmological constant the value y_ for the minimum must be negative. 
Again, the limiting case (y_ = 0) corresponds to fi = (J>~{\) as can be seen in Figures 14 
and 12. In any case the necessary and sufficient condition for the cosmological constant to 
be real is y_(oo) < or j/jv(oo) < h. Therefore we have an excluded region of parameters 
below fi = The other curve \i = /i + (A) does not affect the qualitative behaviour of 
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Figure 14: A(r) < V?* £ [r + ,r^ ) U (r+,00) corresponding to M^ 2 \ In this case, as in Fig. 12, 
there is no way of continuously connect the real cosmological constant to the horizon. 



the relevant solution since it is just indicating the appearance of two new real cosmological 
constants. Thus, in most of the space of parameters the existence of two extra solutions 
does not qualitatively affect the solution with horizon, except in the excluded region where 
the cosmological constant connected to the horizon is not real. 
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Notice that the well behaved solution is regular when crossing the curve fi = ^ + (A); 
the would be symmetry enhancement affects the other two cosmological constants. In other 
words, the two cosmological constants that agree over this curve are not those connected 
to the horizon, and then the theory has propagating linear perturbations when expanded 
about our vacuum. Symmetry enhancement for this vacuum occurs at fj, = yT (A) and this 
is the boundary of the excluded region. 9 

The usual parameterization of the black hole solution [31, 39] shall make manifest the 
different properties portrayed in the previous section. In fact, the usual parametrization of 
the three, in general complex, solutions for the function / in third order Lovelock gravity 
is 



1 + an ( J(r) + TOM) V3 + a, ( J(r) - v/0(r)) ^ 



where 



with 



J(r) 2 +T 3 



J(r) 



r 



A 2 



3 fJL 

2A2 



1 



3 ii 2 
2A3 



3{i 2 
2A3 



VN(r) 



4 A 6 



/r 



(4.16) 
(4.17) 

(4.18) 
(4.19) 



on are the three cubic roots of unity, ao = 1 , a± = — ; and the bar indicates complex 
conjugation. Each of these solutions is associated with one possible value of the cosmologi- 
cal constant, and so, fixing the value of the cosmological constant fixes the function. Notice 
that / is directly related to our previous parameterization as 



L 2 



(4.20) 



where x is the relevant branch (a, /3 or 7) in each region: 7 for /i < and, whenever \x > 0, 
a for A < and f3 for A > 0. 



5. Causality constraints 

We will now explore constraints arising from the demand that boundary signals propagate 
with velocities lower than the speed of light. This strategy has been originally pursued 
in the case of Gauss-Bonnet [5, 6] and extended to all possible helicities [7, 8] and higher 
dimensions [9, 10]. We will follow, in particular, the conventions used in [10] and we refer 
the reader to that paper for futher details. 

9 If we want to study perturbations about the vacuum at one of these points, from the equations of 
motion we see that 

e afl ... fd _ 1 (V l/2 - A e /l/2 ) A (V 3/4 - A e /3/4 ) A 8 (R hh - A e /5/e ) A e ^'^-i = , 
is a trivial identity and there are no quadratic perturbations. 
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5.1 Black hole perturbations 

We shall consider perturbations of the metric around the black hole solution obtained 
in the previous section along a given direction parallel to the boundary (say, and 
propagating towards the interior of the geometry. Using the direction x d ~ x as an axis of 
symmetry, we can classify the perturbations in helicity representations of the rotation group 
around it. It is convenient to analyze each case separately. The linear order contribution 
to the equations of motion (2.4) can be written as 



S£ a = e a/l ... /d _ 1 °k [kSR flh A Rf*"fv> A e /MH-i~/«j-i 
k=0 

+ (d-2k- l)R fl - f2k Ae /2fe+1 '" /d - 2 Ade 1 "- 1 =0. (5.1) 

The first thing we have to realize in order to carry out this calculation is that the relevant 
contributions to order co 2 , q 2 and u q, come from derivatives along the directions e° and 
e d_1 . In the simplest case, for helicity two perturbations, 10 these contributions have only 
an effect on the expressions of Su 02 , 8u 03 , 5^ d ~^ 2 and 5u^ 1 ^ . Since we are at the 
linearized level, we conclude that the only non-trivial contributions to order u 2 , q 2 , ojq 
come from their exterior derivative (the second term in (5.1) is also irrelevant), 



co 2 , n o uqL 

7T- 6 e he -\ —. = 

2N 2 # f V 2rN#VJ 



5R 02 « d(5co 02 ) = - 7T f^ cp e° A e 2 + n " T „ <f> ^ A e ' , ( 5 - 3 ) 



5R^ 2 « d(6u^ 2 ) = e^ 1 A e 2 - ct> e° Ae 2 , (5.4) 

2r z 2rN#y/j 

where the symbol ~ refers to those contributions relevant to compute the propagation 
speed of a boundary perturbation. There are analog expressions (with the opposite sign) 
for the components with a leg along the e 3 direction. Notice that the cosmological constant 
term is irrelevant for our purposes. 

Recalling from (3.5) that the curvature 2-form of the black-hole space-time is propor- 
tional to e a A e b , it is immediate to see that the only non- vanishing contributions are those 
with d(5uj ab ) behaving equally, 



(and, of course, an analogous expression for d(5u 03 ) and d(5u^ d 1 * >3 )). Some of the equa- 
tions of motion are trivially satisfied. Indeed, notice that 5£ a ~ trivially, unless a = 2 or 

10 We can choose the helicity two perturbation simply as hzsit, r, x 11-1 ) dx 2 dx i . Since we will consider 
small perturbations, h,2z{t,r,x d ~ 1 ) = e<f)(t,r,x d ~ 1 ) with e <C 1. This amounts to e a = e a + e5e a [10], 

e° = N& \/~f dt , e 1 = —5= dr , e F = — dx F , F=A...d-i , 

^7 L 

e = T - (l + | ^) {dx 2 + dx 3 ) , e 3 = ~ (l - | c^j (dx 2 - dx 3 ) . (5.2) 
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3. This is due to cancellations of contributions coming from d(5io 02 ) and d(5uj 03 ) (similarly, 
d(6u}^ d ~^ 2 ) and d(&i/ d-1 ) 3 )). Moreover, by symmetry, the two non-trivial equations just 
differ by a global sign. We must then focus on a single component of (5.1), say, 683 = 0. 
The contribution to first order from each Lovelock term can be written as 



SS 3 



[^1 



e 30 2/ 3 .../ !j _ 1 Yl k Ck d ( duj02 ) A R h "' hl A e hk + 1 - fd - 1 



k=l 



[^1 



'2fe + l---/d-l 



k=l 



. (5.6) 



To proceed, the only thing to worry about is where the and 1 indices are, since depending 
on them the curvature 2-form components change (3.5). Notice that the first (second) line 
gives the to 2 (q 2 ) contribution. The former, for instance, can be nicely rewritten in terms 
of g(r), 

2 K 

SS 3 , U ~ -4r-r <P E k c * Afc "' i^ 1 )' + (d- 3)/" 1 ) (d - 4)! , (5.7) 



while the latter reads 



Nlf 

# J k=0 



k=0 



+ 2(d - 3)r(g k ~ 1 )' + (d - 3)(d - 4)/" 1 ) (d - 5)! . 



(5.8) 



In both expressions ~ signals the ellipsis of an irrelevant product of (d— 1) vierbeins. Here 
we assume that c k >K = 0. It is convenient to define the following functionals, C d [f,r], 
involving up to /cth-order derivatives of / (or g): 



K 



4° ) b,r]=^fcc fe A fe - 1 /- 1 , (5.9) 
k=i 

K 

C^[g,r] = Y,kc k A k - 1 {r(g k ~ 1 y + (d-3)g k - 1 ) , (5.10) 
fc=i 

K 

Cf [9, r] = ^kc k A k ~' (r 2 {g k - 1 )" + 2(d - 3)r(/^)' + (d - 3)(d - 4)/" 1 ) .(5.11) 



fc=i 



Then, the speed of the helicity two graviton can be written as (see [10] for details) 

9 Cd ] [g,r] 



4 



(d-4) c W 



(5.12) 



d 19, r\ 



and expanding about the boundary 



g(r) ^l-k 



„d-l 



,d-l 



(5.13) 
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we get 



i 



1 + 



2(fc— 1) 
(d-3)(d-4) 



„d-l 



(5.14) 



In [5, 6] it has been argued that the existence of a maximum in the radial dependence 
of this speed where it exceeds unity indicates the possibility of bouncing high-momentum 
gravitons leading to causality violation in the boundary theory. Their arguments extend 
smoothly to our case, as discussed in [10] . Now, from the polynomial equation for g we get 
(3.19) and thus, 



4 



i + 



,fc-i 



i + 



2(fc-l)(d-l) 
(d-3)(d-4) 



L 2 A r d ~ l 



£f=i^ fe A^ 



(5.15) 



It is a bit more involved but the same can be done for the other helicities (see [10] for 
the details). For the helicity one case we have to solve the equations of motion for two 
components, but one of them always vanishes. The remaining equation yields directly the 
speed of the graviton as in the previous case 



4 



9 C^[g,r] 



(5.16) 



Expanding again near the boundary 

1 rt'ZtikckA^ 



1 + 



L 2 A r d ~ x 



££=1*<*A*-1 



(fc-i)(rf-i) 

(d-3) 
2 



(5.17) 



The more complicated case is the helicity zero one, where we have to find the four com- 
ponents of the perturbations. One of them [<j)\ see [10]) is set to zero by the equations of 
motion, while other two can be written in terms of only one degree of freedom 



(d-3) £ 



(5.18) 



When we substitute this into the equations of motion, only one of them remains linearly 
independent and gives the speed of the helicity zero graviton, 



9 



' 2 C£Vj _ C^[g,r) 



(5.19) 



that near the boundary yields 



1 + 



1 rt'EtikckA^ 



2(fc-l)(d-l) 
(d-3) 



L 2 Ar d-l 



£f=i^ fe A*-i)' 



(5.20) 
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We shall analyze now the would be restrictions imposed on the Lovelock couplings to avoid 
superluminal propagation of signals at the boundary. The above expressions are valid for 
arbitrary higher order Lovelock theory and higher dimensional space-times. As we saw 
earlier, though, there could be additional restrictions coming from the stability of the AdS 
vacuum solution or, further, the existence of a black hole with a well defined horizon in such 
vacuum. Instead of analyzing these expressions right away, let us discuss an alternative 
computation, first introduced in [8] , given by the scattering of gravitons and shock waves. 



5.2 Shock waves 

The shock wave calculation, originally introduced by Hofman in the case of Gauss-Bonnet 
theory in 5d [8] , and extended by us to arbitrary higher dimensional Gauss-Bonnet gravity 
[10], can also be generalized to general d space-time dimensions in the case of higher order 
Lovelock dynamics. The relevant shock wave solution is w = a z d ~ 3 . The proceedure is 
almost the same as in the previous section, just a bit more complicated since the symmetry 
of the background is lower than in the black hole solution. As before, since we are only 
interested in the high momentum limit, we keep only contributions of the sort d 2 (j), d u d v (j) 
and d 2 (p. These contributions come again from the exterior derivative of the perturbation 
of the spin connection. In the helicity two case 



02\ 



d(5uj 



L 2 z 2 
N l 

L 2 z 2 
N l 



d 2 (f) e 1 A e 2 + ( d u d v <(> + af(u)L 2 z a - 1 d$<j) ) e u A e 



2 d-l o2 , \ „0 A „2 



d u d v <f> + af(u)L 2 z d ~ l d 2 ^) e 1 A e 2 + (• • • ) e u A e 2 , (5.22) 



(5.21) 



the ellipsis being used in the second expression since the corresponding term does not 
contribute to the equations of motion. The components with index 3 instead of 2 are the 
only remaining ones non- vanishing, and they are obtained changing — > — <j). The other 
thing we need is the curvature 2-form of the background metric, that can be written as 



R ab = A(e a A e b + X ab ) , 



(5.23) 



where A = — L -^ N 2 an d X ab is an antisymmetric 2-form accounting for the contribution of 

# 

the shock wave 



X 



la 



{d- l)af(u)L 2 z d - 1 e° Ae" , o^O.d-l, 



X Ud-i) =_[(«*_ 2 ) 2 - 1] a f(u) L 2 z d ~ x e° A 



4-1 
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Now, the relevant equation is, as before, given by a single component of (5.1), 



A' 



5£ 3 = ^2kc k e 3fl ... fd _ 1 d(5uj flh ) AR h - f2k A e f2k+v " fd - 



k=l 
A 



2]Tfcc fc [e 3 02f 3 -f d ^d(Suj 02 ) + e3 i2/ 3 .../ £i -i^ 12 )] A i2 /s "' /2fc A e Afc + 1 "'^- 1 
fe=i 

4L 2 z 2 ' 



^fccfcA*" 1 [(d-3)!(9 u ^ + a/( U )L 2 Z d - 1 9>j -(fc-l)x 
# fc=i 

x(d-5)!((d-5)(d-l)-[(d-2) 2 -l])a/(u)L 2 z^ 1 a 2 (/.j , (5.24) 

since the diagonal part of d{5ui ab ) contributes to all the diagonal parts of R ab and the off- 
diagonal part of d{5uj ab ) contributes to the off-diagonal part of one of the R ab . Collecting 
terms type d u d v cp and d%4> and simplifying constant factors (not depending on the degree 
of the Lovelock term, k), 5£ 3 = implies 



g k c k A*" 1 \dM + (l + 2 y_l][*_q ) V" 1 ^ = . (5.25) 

Causality problems appear when the coefficient of d^cp becomes negative. For the corre- 
sponding values of the couplings, a graviton that is emitted from the boundary comes back 
and lands outside its own light cone. We shall impose, thus, 



(d-3)(<f-4) 



> (5.26) 



This condition coincides with the one obtained in the helicity two perturbation analysis of 
the AdS black hole solution (5.15) provided R > 0, or equivalently , 



A 

J^CfcfcA^X), (5.27) 
k=i 

which was shown to be the case in (3.18). Needless to say, the other two helicities match 
these expectations as well. Then we have shown how the calculation involving bouncing 
gravitons in a black hole background and those scattering shock waves lead exactly to the 
same result for completely general Lovelock gravities in arbitrary dimensions. The general 
constraints can be written as 

helicity 2: r'[A] + -^- A1»[A] > , 

helicity 1 : T'[A] - fc-^AT"[A] > , (5.28) 

(d-3) 

helicity : T'[A] - ^£^AT"[A] > , 
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where T[A] = is the polynomial equation for the cosmological constant. These constraints 
are also valid for a / (different horizon topologies) as in their derivation everything 
depends on / through the universal function g. This independence of the topology of the 
black hole horizon is connected to the alternative calculation using shock waves. 

We would like to show, for definiteness, the results corresponding to third order Love- 
lock theory. In this case, the analysis of the previously discussed causality constraints 
reduce to studying the sign of the following set of polynomials (x = L 2 h))\ 



M 2 
M 

A/n 



2(d 2 -5d+10h d 2 -3d + 8 2 



(d-3)(d-4) 
4 



(d-3)(d-4)' 



1 



(d-3) 
2(d+l) 



Xx 



d+1 2 



3d-l 2 

AX — — —fl X 



(5.29) 
(5.30) 
(5.31) 



(d-3)'"" (d-3)' 

where N stands for the numerator of the first correction to the speed at which perturbations 
propagate at the boundary and the subscript indicates the corresponding helicity. The 




Figure 15: In d — 7, the blue shaded region correspond to negative values of A2 for the branch 
of solutions with a horizon, and therefore causality violation occurs there. The red and blue thick 
curves correspond to fl = ^ [189A + 4 (—8 ± (2 — 9A) \/16 — 45A)] . The blue one provides the 
restriction for the relevant solution branch. The thin blue and red lines are nothing but the singular 
locus A = 0. We can observe that this region includes the Gauss-Bonnet case (Aqb = 3/16, fl = 0) 
(depicted as an orange point) [9], If fi = A 2 , contrary to what is stated in [39], there is a forbidden 
region limited by the left green point, A > 19/81 (the right one being the maximally symmetric 
(Chcrn-Simons) point, A = 1/3). 



constant x, as we discussed earlier, is one of the solutions of 

A* 
3 



l + x + Xx 2 + ^x 3 = 0. 



(5.32) 
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To proceed, we have to find the simultaneous solutions for each of the polynomials and the 
previous equation, that is where each of the polynomials can change sign. The simplest 
way of doing this is to find the solutions for each of the polynomials and substitute them 
into the equation. This procedure gives several curves [i = fi(X) for each polynomial. The 
most clear way that we find to present these results is by presenting the following figures 
which illustrate the situation. The analytic expression of the curves can be easily obtained 
and, indeed, are discussed in the captions. We will discuss in the text their main qualitative 
features. 

In Figure 15 we display both the region of the coupling space excluded by the arguments 
of the preceding section and (in blue) the sector of the (A, //) plane that leads to causality 
violation in the 7d Lovelock theory. Projected in the A axis we reobtain the upper bound 
corresponding to Gauss-Bonnet gravity and obtained in [9, 10]. Even if, from the point 
of view of AdS/CFT, the trustable region is close to the origin, i.e., for small values of 
A and fi, it is amusing to see that the forbidden region has a structure that seems to be 
meaningful, even if this is not fully justified. Notice that, had we restricted to third order 
Lovelock theory with fi = A 2 , we would have obtained a higher upper bound for A (this 
corrects a misleading statement in [39]). 

We represented the situation corresponding to helicity one gravitons colliding a shock 
wave or, alternatively, helicity one perturbations of the AdS black hole background in 
Figure 16. We see that it gives a complementary restriction in the coupling space. Again, 



Figure 16: The helicity 1 polynomial restrict the parameters to lie inside the region between 
the blue curve and the excluded region. The shaded region correspond again to negative values 
of Mi for the relevant branch and thus to causality violation. The equation of the curve is fi = 
[-18A + 49 ± (7 + 6A)\/49 - 120A] . For Gauss-Bonnet we found, A > -2. 

it intersects the A axis at the point —2, as obtained in [10]. We will not extend much in 
the discussion of this curve for reasons that will be clear immediately. 
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Let us finally explore the causality violation induced by helicity zero gravitons. A very 
important feature is the fact that this curve imposes more severe restrictions than those 
forced by the helicity one case. This can be easily seen to be a general statement arising 
from the constraints (5.29). Indeed, the curve fully overcomes the helicity one constraints 
in such a way that they are finally irrelevant. This has been seen earlier in Gauss-Bonnet 
[8, 10]. Going to the conformal collider physics setup, this result seems to be related to the 
vanishing of £4 (see the Appendix). The intersection with the A axis reproduces the result 




Figure 17 : As in the helicity 1 case, the helicity constraint restrict the parameters to lie in a strip 
between the blue line and the excluded region, but it is always more stringent. For Gauss-Bonnet, 
A > The equation of the curve is M = nk [ 315A + 128 ± 4 ( 4 + 15AV64- 165AJ . 

for Gauss-Bonnet [9, 10] (see Figure 17). 

Now, the three curves should be put altogether to determine the region in the coupling 
space that is allowed by all helicities gravitons. As mentioned, the helicity one curve ends up 
being irrelevant. The coupling space of the third order Lovelock theory becomes, as we see 
in Figure 18, quite restricted. Several questions immediately arise and we will try to pose 
some of them and answer quite a few in the concluding section. It should be pointed out 
that there is a subregion of the allowed sector with A greater than 1/3 which immediately 
leads to negative values of 77/s [25]. This result is unphysical and clearly deserves further 
study. 

To finish this section, it is tempting to use our expressions for third order Lovelock 
theory which are valid for arbitrary higher dimensional space-times. It is not difficult to 
compute the allowed regions for different dimensionalities. This is represented in Figure 
19. It should be pointed out that the subregion which, being part of the allowed sector, 
leads to negative values of the shear viscosity disappears for lid (and higher). Indeed, the 
maximum allowed values of A are, respectively, 1369/3519, 147/377 and 2209/5655, for 8d, 
9d and lOd, which in all the cases exceed the critical value leading to a vanishing shear 
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Figure 18: As the helicity constraint is always more stringent than the helicity l's, the allowed 
region of parameters is contained between the curves corresponding to helicity and 2 (red and 
blue lines respectively). Notice that the maximum value for A is raised from the GB case to the 
value 64/165 (yellow point). Strikingly enough this gives a negative value for the shear viscosity to 
entropy density ratio. 

viscosity (respectively, 5/14, 3/8 and 7/18). As stated, lid is the first dimensionality at 
which the maximum value for A, 169/432, is lower than the critical one, 2/5. Whether 
this is a numerical coincidence or it is pointing out that something special happens to 
(super)gravity in lid is presently unclear. 

6. Concluding remarks 

Let us start by summarizing the contents of this article. We have used the AdS/CFT 
framework to scrutinize higher order Lovelock theory in arbitrary dimensions with the 
focus on possible causality violation. We have formally addressed the more general case 
and explicitely worked out the third order Lovelock gravity in arbitrary dimensions. Results 
for higher order theories are contained in our formulas though some extra work is needed 
to extract them explicitly. The third order case has taught us that it is a subtle issue 
to determine which is the branch of solutions connected to Einstein-Hilbert and, thus, 
presumably stable. From the point of view of a black hole background, the relevant question 
is whether one can find an asymptotically AdS black hole with a well-defined horizon. This 
is a delicate problem that can be amusingly casted in terms of a purely algebraic setup. 

We do this computation in two different ways and we end up showing that they are 
bound to give the same result. We have first computed all polarization linear perturbations 
of the black hole AdS solution in Lovelock theory, and subsequently studied the collision 
of gravitons and shock waves in AdS in the framework of this theory. We found a region 
in the coupling space where the theory does not violate causality (at least in a way that 
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0.6 r 




Figure 19: Allowed region for the GB and cubic Lovelock couplings for arbitrary dimension. The 
inner region is the one depicted in Figure 18 for d = 7. Towards the exterior, we depicted with a 
decreasing blue shadow successively, those for d = 8, 9, 10, 11, 20, 50 and oo. Blue points denote the 
boundaries of d = 5 and 6, where the cubic Lovelock term identically vanishes. The orange points 
are the asymptotic values for GB gravity found in [10]. Notice that one of the branches that bound 
the allowed region tend to the singular locus A = when d — > oo. 

can be detected by this kind of computation). This generalizes the intervals found for A 
in the case of Gauss-Bonnet theories [7, 8, 9, 10]. Contrary to those intervals (that are 
now just given by the intersection of the allowed sector and the A axis in this paper), 
the region is unbounded. Arbitrary large (simultaneous) negative values for A and \x are 
allowed. The meaning of this result is unclear since this happens in a region of couplings 
where the computations are not trustable. However, the regular pattern that we found as 
a function of the space-time dimensionality is too nice to be discarded so easily. We have 
seen that the allowed region increases with space-time dimensionality, as it happened with 
the Gauss-Bonnet coupling [10]. In the infinite dimensional limit, one of the boundaries of 
the allowed region (roughly speaking, the one giving the upper bound; rigorously speaking, 
the one due to helicity two gravitons), asymptotically approaches the discriminant of the 
Lovelock theory. This is expected since the region beyond that curve is excluded. The 
other curve, arising from scalar gravitons, has a less natural behavior that calls for a 
deeper explanation. There is no restriction for this curve to be asymptotically larger, but 
it approaches a finite curve, that of course crosses the A axis at A = —3/4, in accordance 
with [10]. Thus, moderate (non- negative) values of A and [i are forbidden even for would 
be infinite dimensional theories. Whether this is an artifact of our approach or there is a 
deep reason for this behavior, is unclear to us. 

The papers that originally dealt with causality violation did it under the spell of 
searching possible violations of the KSS bound [5, 6]. They found that the addition of 
a Gauss-Bonnet term in the gravity Lagrangian corrects the universal value of the shear 
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Figure 20: We have shadowed the region of the allowed space of couplings leading to unphysical 
negative values of rj/s. It can be seen that the maximum A increases with the space-time dimen- 
sionality, but the critical value leading to a vanishing r\j s does it faster. From lid on, this problem 
disappears. 

viscosity to entropy density ratio as 



The addition of further higher order Lovelock terms does not manifestly contribute to this 
quotient [24]. The possible addition of a quadratic curvature correction with positive A is 
enough to argue that the KSS bound is not universally valid. The appearance of an upper 
bound for A due to causality conservation naively seems to lower the KSS bound down to 
a new one. Since (6.1) is not affected by higher order Lovelock terms, this has led some 
authors to speculate about the nature of this newly found lower bound. However, it is not 
clear if the new bound should be taken seriously: when A approaches the maximal allowed 
value, the central charges of the dual CFT explore a regime (roughly (a — c)/c ~ 0(1)) 
where the gravity description is a priori untrustable. Furthermore, let us point out that 
the fact that higher order Lovelock couplings do not enter (6.1) does not mean that they 
are irrelevant for this problem. Indeed, we have seen that even an infinitesimal positive 
value of [i would raise the upper bound for A thus affecting the putative lower bound for 
rj/s. Even more weird, there are values of fi leading to such high values of A that -q/s 
becomes negative. There should be some other mechanism or argument to disregard this 
unphysical behavior. 

In previous articles on the Gauss-Bonnet case it was raised the issue about the string 
theory origin of quadratic curvature corrections. It seems clear after [15, 16, 17] that these 
terms may possibly arise from D-brane probes or A^^i singularities in M-theory. The 
situation with the cubic curvature corrections is more delicate in this respect. Besides, 
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the very fact that the results generalize smoothly both to higher order Lovelock theory 
and to higher space-time dimensions, seems to suggest that these computations may not 
necessarily rely in the framework of string theory. 

Our approach lacks from the computation in the field theory side adapted to the 
framework of conformal collider physics [2]. This makes us difficult to explicitly relate the 
restrictions imposed by causality, obtained in this paper, with those presumably emerging 
from positive energy considerations. We would need to know, for instance, how the higher 
order Lovelock coefficients affect the values of the dual CFT central charges. In any case, we 
can make some qualitative statements. Our results for arbitrary dimensions are compatible 
with the condition £4 = 0. This is the expected value of £4 for a supersymmetric CFT. 
Cubic terms are expected to accomplish supersymmetry breaking. However, this is not the 
case for the Lovelock combination. As for t2, the result of this paper is unexpected. It is 
well-known that the third order Lovelock term endowes a 3-graviton vertex structure that 
does not contribute to the 3-point function [22]. However, this result is strictly valid for 
flat space-times while in AdS, as we see, it is not true. The introduction of the third order 
Lovelock Lagrangian enters non-trivially in the discussion of causality violation. It would 
be interesting to further clarify how this term contributes to £2 by studying directly the 
3-graviton vertex in AdS space-times [40]. 
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A. Conformal collider physics 

We summarize in this appendix a few formulas and explanations that we need in the bulk 
of the paper regarding the so-called conformal collider setup [2] as applied to conjectural 
higher dimensional CFTs. We would like to measure the total energy flux per unit angle 
deposited in calorimeters distributed around the collision region , 



the vector h pointing towards the actual direction of measure. The expectation value of 
the energy on a state created by the stress-energy tensor, O = eijTij, thus, is given in 
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terms of 2- and 3-point correlators of [2], 

(0|Otg(n)O|Q) 

= (o|oto| ) • (A - 2) 

Using the fact that €ij is a symmetric and traceless polarization tensor with purely spatial 
indices, this is fixed up to coefficients t 2 and £4, 



(£{n)) e T = 



1 I t 2 ( ^ ^ ^ Hj — I +t 4 ' l^'™*"^ 2 2 



d-2y ^ e*^ d(d-2)y 

■ (A.3) 

the d dependence in this expression being given by the normalization of the integrals over 
d — 3 sphere whose volume is 0^-3- By demanding positivity of the deposited energy 
irrespective of the calorimeter angular position, a conformal collider gedanken experiment 
in a higher dimensional CFT would yield positive energy bounds from (A.3): 

tensor: 1 - t 2 - U > , (A.4) 

vector: (l - ^ t 2 - ^-^ U ) + \t 2 > , (A.5) 

scalar : (l - ^ t 2 - t 4 ) + ^| (t 2 + U) > . (A.6) 

The three expressions come from the splitting of into tensor, vector and scalar compo- 
nents with respect to rotations in the hyperplane perpendicular to h. These constraints 
restrict the possible values of t 2 and t± for any CFT to lie inside a triangle whose sides are 
given by (A.4), (A.5) and (A.6). The vertices of the triangle are (— ^-sd+i ' 331)1 (0) ^2 ^ ) 
and (d,—d). It is a straightforward consequence to see that the helicity one contribution 
is not restrictive for £4 < -r^y. In particular, of course, it is so for £4 = 0. 
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